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I.  INTRODUCTION  AND  SUMMARY 


1.1  INTRODUCTION 

The  objective  of  the  Systems,  Science  and  Software  (S3) 
research  program  is  to  extend  our  present  understanding  of  the 
generation,  propagation,  detection  and  analysis  of  seismic 
waves  by  both  underground  explosions  and  earthquake  sources. 

The  goal  is  to  improve  the  United  States’  ability  to  monitor 
compliance  with  treaties  limiting  underground  nuclear  explo¬ 
sions.  This  report  covers  work  performed  at  S3  during  the 
twelve  month  contract  period  between  1  October  1979  and 
1  October  1980,  under  Contract  Number  F08606-79-C-0008  (VSC 
Project  VT/0712) . 

The  report  includes  six  distinct  technical  sections, 
which  are  briefly  summarized  below,  together  with  the  refer¬ 
ences  to  the  appropriate  task  in  the  Work  Statement  of  the 
contract.  Much  of  this  work  has  been  described  previously 
in  several  S3  topical  and  technical  reports.  In  such  cases, 
we  have  provided  an  executive  summary  of  the  work  performed 
during  the  contract  period,  including  a  digest  of  the  principal 
conclusions  and  the  relevant  figures. 


1.2  SUMMARY  OF  SECTION  II:  SEISMIC  SOURCE  CALCULATIONS 
This  section  contains  three  subsections. 

•  Development  of  a  constitutive  model  for  salt 
based  on  an  analysis  of  free-field  ground 
motion  data  for  nuclear  and  high  explosive 
events  in  salt.  The  events  included  GNOME, 
SALMON,  STERLING  and  COWBOY.  Because  of  the 
large  body  of  digitized  data  available  for 
SALMON,  this  event  was  chosen  for  normaliza¬ 
tion  purposes.  The  purpose  was  to  explain 
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simultaneously  all  observations  with  a  single 
constitutive  model  for  salt  so  as  to  permit 
development  of  realistic  decoupling  scenarios. 
Good  agreement  with  data  from  STERLING  and 
GNCME  was  produced.  For  COWBOY,  predicted 
and  observed  peak  velocities  were  in  agree¬ 
ment,  but  predicted  peak  displacements  were 
larger  than  those  observed  by  a  factor  of 
two.  This  work  was  done  by  Dr.  J.  Theodore 
Cherry  and  Dr.  Norton  Rimer. 


•  Analysis  of  the  dynamical  aspects  of  seismic 
sources  which  yield  different  signatures  for 
earthquakes  and  explosions  was  conducted  through 
numerical  simulations  of  earthquake  sources. 

The  simulations  were  performed  on  the  ILLIAC  IV 
computer,  using  a  three-dimensional  finite 
difference  code.  Inelastic  response  along  the 
fault  zone  was  modeled  through  nonlinear  boundary 
conditions  on  the  fault  plane.  As  a  result  of 
this  work,  it  appears  necessary  to  include  more 
complex  rupture  physics  in  the  modeling  procedure, 
and,  in  particular,  a  better  description  of  the 
cessation  of  rupture  is  required  to  explain  the 
discrepancies  between  modeled  and  observed 
stopping  phases.  Stress  irregularities  have 
been  shown  to  have  a  significant  influence  on 
the  spectral-  and  time-domain  characteristics 
of  the  teleseismic  signals  from  earthquakes;  they 
tend  to  increase  the  separation  of  earthquake  and 
explosion  populations  for  the  m^-Mg  discriminant. 
This  work  was  performed  by  Dr.  Steven  M.  Day. 
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•  A  simulation  study  for  a  (hypothetical)  fully 
decoupled  5.3  KT  explosion  at  the  SALMON  site 
was  conducted  by  scaling  regional  phases  re¬ 
corded  in  the  eastern  United  States  for  SALMON. 
Such  an  event  could  probably  be  detected  at 
epicentral  distances  up  to  ^  650  km,  but  not 
at  distances  greater  than  ^  1000  km.  For  low- 
noise  conditions,  the  Lg  phase  would  probably 
be  detected  at  distances  up  to  'v  1500  km,  but 
not  the  initial  P  phase.  This  study  was  per¬ 
formed  by  Mr.  John  R.  Murphy  and  Dr.  T.  Joseph 
Bennett. 


1.3  SUMMARY  OF  SECTION  III:  REGIONAL  SEISMOLOGY 

This  section  pertains  to  a  study  of  far-field  seismic 
waves  associated  with  a  series  of  two-dimensional  calculations 
of  buried  explosions.  The  purpose  was  to  determine  the  impor¬ 
tance  of  nonlinear  free  surface  interactions,  including  spalla¬ 
tion,  and  the  depth  dependence  of  overburden  pressure  and 
rheology.  These  effects  result  in  higher  order  source  prop¬ 
erties,  which  constitute  a  departure  from  the  classical 
spherically  symmetric  point  source  approximation.  This  work 
has  been  the  subject  of  a  detailed  report  by  Drs.  Thomas  C. 
Bache,  Terrance  G.  Barker,  Norton  Rimer  and  J.  Theodore  Cherry 
(1980)  and  we  only  include  the  executive  summary  for  complete¬ 
ness  . 


1.4  SUMMARY  OF  SECTION  IV:  AUTOMATIC  SIGNAL  PROCESSING 
This  section  includes  two  subsections. 

•  A  review  of  the  MARS  (Multiple  Arrival 
Recognition  System)  detector.  The  MARS- 
based  detection  algorithm  was  optimized 
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and  an  efficient  code  was  tested  on  an  off-line 
basis  on  data  supplied  by  VELA  Seismological 
Center  (VSC) .  This  subsection  was  written  by 
Dr.  William  E.  Farrell,  Mr.  Robert  C.  Goff  and 
Mr.  James  Wang. 

•  A  review  of  state-of-the-art  detection  algorithms, 
by  Dr.  Jonathon  Berger.  This  section  identifies 
and  evaluates  all  pertinent  automatic  detection 
algorithms,  and  provides  a  classification  by 
equation,  mathematical  derivation,  and  usage 
history.  It  includes  algorithms  currently  under 
development  and  scheduled  to  be  validated  by 
October  1981.  In  addition,  recommendations  for 
further  research  efforts  are  made. 


1.5  SUMMARY  OF  SECTION  V:  EVALUATION  OF  SUITES  OF  EVENTS 

This  section  summarizes  work  by  Dr.  Thomas  C.  Bache 
and  Mr.  Boris  Shkoller  which  is  described  more  fully  in  a 
semiannual  technical  report  by  Bache,  et  al.  (1980) .  It 
pertains  to  the  automatic  determination  of  body  wave  magni¬ 
tudes  using  the  MARS  processor,  and  the  testing  of  a  similar 
automatic  surface  wave  magnitude  estimator.  Particular  atten¬ 
tion  is  paid  to  the  effect  of  seismic  noise,  which  are  evalu¬ 
ated  by  superposition  of  a  signal  with  known  magnitude  and 
seismic  noises  of  varying  levels.  The  results  appear  to  be 
robust  with  respect  to  noise  even  in  the  case  of  high  noise 
level. 

Applications  to  actual  records  from  Eurasian  events  at 
RKON  have  also  been  performed  and  have  been  described  in  the 
aforementioned  report. 
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1.6  SUMMARY  OF  SECTION  VI:  GROUND  MOTION  ANALYSIS 

This  work,  by  Mr.  John  R.  Murphy  and  Dr.  T.  Joseph 
Bennett  pertains  to  the  compilation  and  evaluation  of  near¬ 
source  data  from  underground  explosions.  In  particular,  an 
assessment  of  the  data  set  in  terms  of  the  seismic  source 
function  for  different  emplacement  media  has  been  conducted. 
The  expanded  data  set  includes  previously  unpublished  data 
from  the  Sandia  Corporation  archives.  A  detailed  report  has 
been  provided  by  Murphy  and  Bennett  (1980). 

1.7  SUMMARY  OF  SECTION  VII:  ACOUSTIC  WAVES  FROM 

EXPLOSIONS  IN  THE  OCEAN 

The  results  for  this  section  have  been  reported  by 
Bache,  et  al.  (1980) .  They  specifically  include  estimates 
of  the  gross  spectral  characteristics  and  signal  duration  of 
acoustic  waves  generated  by  a  nuclear  explosion  near  the 
ocean  surface.  A  nominal  yield  of  1  KT  and  nominal  range  of 
6600  .cm  are  assumed. 
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II.  SEISMIC  SOURCE  CALCULATIONS 


2.1  INTRODUCTION 

Investigations  of  seismic  sources  by  theoretical  and 
numerical  calculations  were  covered  by  three  tasks: 

One-  and  two-dimensional  explosion  calculations 
(Task  4.1.1}.  (See  also  Section  II,  Task  4.2.) 

Earthquake  simulation  studies  using  the  ILLIAC  IV 
(Task  4.1.2) . 

Decoupling-evasion  research  (Task  4.1.3). 

Progress  accomplished  under  these  tasks  during  the  contract 
period  is  described  in  the  following  subsections,  where 
specific  references  are  made  to  reports  in  which  greater 
detail  can  be  found. 

2.2  A  CONSTITUTIVE  MODEL  FOR  SALT  BY  J.  T.  CHERRY  AND 

N.  RIMER 

Several  numerical  calculations  of  underground  explosions 
performed  at  S3  have  been  used  principally  in  Task  4.2,  and 
the  results  are  summarized  in  Section  II  of  this  report. 

In  this  subsection,  we  focus  on  calculations  designed 
to  permit  development  of  realistic  decoupling  scenarios  in 
salt.  For  this  purpose,  it  is  necessary  to  understand  free 
field  motion  data  for  actual  explosions  in  this  material. 

We  have  used  ground  motion  data  from  nuclear  and  high  explo¬ 
sive  events  including  GNOME,  SALMON,  STERLING  and  COWBOY. 

Because  the  SALMON  event  provided  a  large  subset  of 
the  data,  and  because  many  of  the  associated  records  are 
available  in  digital  form,  we  chose  this  event  for  the  purpose 
of  normalizing  the  constitutive  model.  This  model  was  then 
applied  to  the  other  events,  in  order  to  determine  how  much 
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of  the  salt  ground  motion  observations  could  be  explained  using 
a  single  constitutive  model. 

This  procedure  produced  good  agreement  with  the  ground 
motion  data  from  STERLING  and  GNOME.  The  predicted  peak  dis¬ 
placements  for  the  tamped  COWBOY  shots  were  approximately  a 
factor  of  two  greater  than  those  observed.  However,  predicted 
and  observed  peak  velocities  for  COWBOY  were  in  close  agree¬ 
ment. 


2.2.1  Normalization  to  SALMON 


Figure  1  shows  the  laboratory  strength  data  for  SALMON, 
GNOME  and  polycrystalline  salt.  The  curve  shown  in  the  figure 
fits  the  data  for  SALMON  and  GNOME  and  was  used  in  the  model 


for  the  ultimate  strength  of  salt  (Y 


lim 


)  • 


The  best  agreement  with  the  SALMON  ground  motion  data 
was  obtained  by  assuming  that  salt  work  hardened  to  its  ulti¬ 
mate  strength.  The  expression  used  to  calculate  the  strength 
(Y)  as  a  function  of  the  inelastic  energy  (E)  deposited  in 
the  material  during  yielding  was 


Y0(l  +  ex  E  -  e2E‘)  <  *Iim 


where  YQ  is  the  initial  (saturated)  strength  and  e^  and  &2  are 
respectively  work  hardening  and  work  softening  material  con¬ 
stants. 

If  is  the  deviatoric  component  of  the  stress  tensor 
and  J2'  is  the  unadjusted  second  deviatoric  invariant,  then 
the  increment  in  inelastic  energy  becomes 


(2) 
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Strength  data  for  salt 


where 


3  J'  >  Y2 

it  .  .  »  the  deviatoric  component  of  the  inelastic 
strain  tensor 

U  =  the  shear  modulus. 


This  model  is  simple  and  easily  implemented  in  a  finite 
difference  stress  wave  code.  The  physical  basis  for  the  work 
hardening  aspect  of  the  model  may  be  an  increase  in  effective 
stress,  and  conversely  a  decrease  in  pore  fluid  pressure, 
caused  by  dilatancy.  If  this  is  the  case,  then  YQ  is  the 
saturated  material  strength.  During  shock  loading  microcracks 
open,  the  pore  fluid  pressure  decreases  and  the  strength  in¬ 
creases.  During  release,  the  cracks  close,  pore  fluid  pres¬ 
sure  increases  and  strength  decreases.  These  modes  of 
material  behavior  are  quantitatively  contained  in  the  model 
with  the  rate  of  hardening  and  softening  determined  by  the 
constants  e^  and 

Figure  2  compares  the  calculated  and  observed  ground 
motion  from  SALMON  at  a  range  of  620  ra  where  the  calculations 
were  performed  with  the  work  hardening  model.  The  emergence 
of  a  precursor  with  a  peak  velocity  of  approximately  0.5  m/s 
is  adequately  modeled  by  assuming  a  saturated  strength  (YQ) 
of  25  bars.  Work  hardening  is  required  in  order  to  match  the 
observed  pulse  shape.  At  a  radius  of  967  m,  the  peak  velocity 
in  the  calculation  is  determined  by  the  precursor.  Elastic 
material  behavior  occurs  beyond  this  radius  in  the  calculation. 

Figure  3  compares  calculated  and  observed  peak  veloc¬ 
ities  and  peak  displacements.  The  results  of  two  calculations 
are  shown  in  this  figure;  one  (dashed  curve)  in  which  Yj  ^ 
equaled  the  maximum  strength  (0.68  kb)  given  in  Figure  1,  and 
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VELOCITY  (CM/SEC) 


M  at  u  u  a.* 


TIME  (SEC) 
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* 


2.  Comparison  between  calculated  (dashed  curve)  and 
observed  (solid  curve)  ground  motion  from  SALMON 
at  a  range  of  620  m. 
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Figure  3.  Comparison  between  calculated  and  observed 

peak  velocities  and  displacements  for  SALMON. 
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the  second  (solid  curve)  in  which  varied  with  stress 

state  (P)  as  shown  in  Figure  1. 

Figure  4  shows  the  RVP  spectra  ( [ f ( )  for  these  two 
calculations.  Also  shown  are  the  spectra  of  Springer,  et  al. 
(1968) ,  which  represents  an  average  of  the  free  field  SALMON 
data,  and  Murphy,  which  was  obtained  by  scaling  the  spectrum 
of  the  GNOME  event  to  the  yield  and  depth  of  SALMON. 

These  results,  involving  comparisons  between  calculated 
and  observed  ground  motions  and  RVP  spectra  indicate  that  the 
constitutive  model  adequately  explains  the  SALMON  data.  We 
now  apply  this  model  to  other  salt  events. 

2.2.2  STERLING 

Simulation  of  the  STERLING  event  was  accomplished  by 
detonating  0.38  KT  in  the  SALMON  cavity.  The  RVP  spectrum  for 
STERLING  was  calculated  and  divided  into  the  calculated  SALMON 
spectra  shown  in  Figure  4  (Y^  .  m  »  Ymay.) .  The  resulting  spec¬ 
tral  ratios  are  shown  in  Figure  5.  In  this  figure  we  also 
show  the  averaged  spectral  ratios  obtained  by  Springer ,  et  al. 
(1968)  from  seismic  stations  common  to  the  two  shots.  Agree¬ 
ment  is  excellent  over  the  frequency  range  analyzed  by  these 
investigators . 

2.2.3  GNOME 

The  GNOME  event  was  simulated  using  the  work  hardening 
model.  Figures  6  and  7  compare  calculated  and  observed  peak 
velocities  and  displacements  for  this  event.  The  calculated 
elastic  radius  was  790  ra. 

Figure  8  compares  the  Mueller-Murphy  estimate  of  the 
GNOME  RVP  spectrum  with  the  calculated  spectrum.  Agreement 
is  surprisingly  good  considering  that  the  Mueller-Murphy 
spectrum  was  based  on  ground  motion  data  at  a  range  of  300  m, 
a  distance  well  within  the  region  of  nonlinear  material  behavior. 


Figure  4.  Comparison  of  calculated  and  observed  RVP  spectra  for  SALMON 


GNOME  563 


RANGE  (m) 


Figure  6.  Comparison  of  calculated  and  observed  peak 
velocities  for  the  GNOME  event. 
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2.2.4  COWBOY 


Figures  9  and  10  compare  calculated  and  observed  peak 
velocities  and  peak  displacements  for  the  tamped  shots  in  the 
COWBOY  series.  Agreement  is  good  for  peak  velocities.  How¬ 
ever,  the  calculated  peak  displacements  are  a  factor  of  two 
larger  than  those  observed. 

There  are  a  number  of  possible  explanations  for  this 
result,  e.g.,  the  COWBOY  salt  was  not  saturated,  the  HE  was 
not  adequately  simulated,  the  peak  displacements  were  not 
free  field,  the  charge  emplacement  was  not  consistent  with  a 
spherically  symmetric  geometry,  the  material  properties 
(strength  and  wave  velocities)  of  the  COWBOY  salt  are  different 
from  GNOME  and  SALMON. 

These  issues  are  worth  pursuing,  but  not  here.  Our 
objective  was  to  develop  a  constitutive  model  for  salt  based 
on  SALMON  ground  motion  data  and  then  to  determine  its  range 
of  applicability  for  other  salt  events,  ’"he  results  of  the 
model  are  consistent  with  data  from  STERLING  and  GNOME.  They 
are  not  consistent  with  the  data  from  the  tamped  COWBOY  series. 

2.3  EARTHQUAKE  SIMULATION  STUDIES  USING  THE  ILLIAC  IV 

BY  S.  M.  DAY 

Differentiation  of  earthquakes  from  underground  nuclear 
explosions  on  the  basis  of  their  seismic  signatures  ultimately 
depends  on  features  related  to  the  details  of  source  dynamics. 
We  have  conducted  a  series  of  numerical  calculations  on  the 
I  ILIAC  IV  computer  designed  to  provide  a  realistic  three- 
dimensional  model  of  the  earthquake  rupture  phenomenon.  These 
are  dynamical  simulations,  which  can  be  compared  in  turn  with 
simple  elastostatic  or  kinematical  solutions  in  order  to  help 
identify  the  effect  of  rupture  complexities  on  the  observable 
seismograms.  This  research  has  been  previously  reported  by 
Bache,  et  al.  (1980a) . 
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Figure  9.  Comparison  of  observed  and  calculated  peak 
velocities  for  the  tamped  COWBOY  events. 


f 


19 


srsTens.  science  and  software 


PEAK  DISPLACEMENT  (mils/lb* *1/3) 


SYSTEMS.  SCIENCE  ANO  SOrTWAHE 


We  are  using  a  three-dimensional  finite  difference 
method  to  model  the  features  of  earthquake  dynamics  believed 
to  be  important  for  predicting  the  teleseismic  signal.  The 
calculations  have  simulated  inelastic  response  of  the  fault 
zone  through  nonlinear  boundary  conditions  on  the  fault  plane. 

The  numerical  method  admits  more  general  nonlinear  material 
behavior,  which  is  being  incorporated  as  required. 

Our  initial  numerical  simulations  were  for  a  simple, 

constant  rupture  velocity,  constant  stress-drop  earthquake 

model.  These  results  show  a  significant  influence  of  fault 

width  on  the  slip  function.  For  points  on  the  fault  which 

are  more  than  a  fault  width  away  from  the  focus,  the  slip 

function  is  spatially  quite  uniform.  The  static  slip  is 
2 

given  by  Aa w/pB  ,  in  agreement  with  estimates  based  on 
elastostatic  solutions.  The  numerical  simulations  predict 
a  rise  time  of  w/2B  and  peak  slip  velocity  of  /2fw/8  Aa/pB. 

In  these  expressions,  w  is  the  fault  width.  Act  the  stress- 
drop,  f  the  maximum  frequency,  8  the  shear  speed,  and  p  the 
density. 

Careful  studies  of  the  broadband  character  of  earthquake 
teleseismic  signatures  indicate  that  the  relative  excitation  of 
long-  and  short-period  waves  is  inconsistent  with  predictions 
from  the  simple  model.  In  addition,  a  few  isolated,  severe 
stopping  phases  dominate  near-field  accelerograms  synthesized 
from  the  simple  model,  in  conflict  with  strong  motion  data. 

More  complex  rupture  physics  is  required  to  explain  the  dis¬ 
crepancies  between  the  model  predictions  and  teleseismic  and 
strong  motion  observations.  In  particular,  a  better  physical 
description  of  the  cessation  of  rupture  growth  appears 
necessary. 

A  displacement-weakening  constitutive  model,  described 
by  Figure  11,  has  been  applied  to  simulate  spontaneous  rupture 
in  the  presence  of  both  uniform  and  heterogeneous  prestress. 
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Figure  11.  Displacement-weakening  constitutive 
model  for  the  fault  zone. 
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In  this  model  of  failure,  the  shear  strength  of  the  fault  zone 
is  assumed  to  be  a  decreasing  function  of  its  deformation; 
energy  is  dissipated  during  rupture  extension  and  shear  stress 
is  bounded  on  the  fault  edges.  When  rupture  initiates  in  a 
limited  region  of  high  stress  (Figure  12a) ,  the  model  leads  to 
spontaneous  propagation,  deceleration  and  stopping  of  faulting. 
Figure  13  shows  the  evolution  of  the  rupture  front  in  this 
case,  as  well  as  the  far-field  P  and  S  wave  displacement  spec¬ 
tra  and  the  corresponding  displacement,  velocity,  and  accelera¬ 
tion  pulses.  In  the  presence  of  multiple  high-stress  regions 
(Figure  12b) ,  each  isolated  stress  concentration  behaves 
roughly  as  an  independent  event.  As  shown  in  Figure  14,  rup¬ 
ture  velocity  becomes  quite  irregular,  and  is  strongly  coupled 
to  peak  velocity  of  slip  on  the  fault  surface.  Figure  15  shows 
the  irregular  rupture  growth,  as  well  as  the  far-field  spectra 
and  time-domain  pulses  which  result.  Predicted  accelerograms 
are  in  better  qualitative  accord  with  strong  motion  data. 

The  model  predicts  significant  consequences  for  the 
spectral  and  time-domain  characteristics  of  the  earthquake 
teleseismic  signal  due  to  the  presence  of  stress  irregularities. 
For  example.  Figure  16  compares  the  short-period  teleseismic 
P-wave  for  the  stress-concentration  model  with  that  for  a  uni¬ 
form  stress  model  (scaled  to  a  source  dimension  equal  to  the 
dimension  of  the  stress  concentration) .  Both  the  period  and 
amplitude  of  the  b  phase  (first  peak-to-trough)  agree  for  the 
two  models,  while  the  stress-concentration  model  has  twice  the 
moment  of  the  uniform  stress  model.  The  two  events  have  been 
superimposed  on  an  m^  versus  Mg  plot  (from  Filson  and  Bungum, 
1972)  to  illustrate  that  the  stress-concentration  phenomenon 
enhances  separation  of  earthquake  and  explosion  populations. 
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STRESS-CONCENTRATION  MODEL  (b)  MULTIPLE  STRESS-CONCENTRATION  MODEL 


Stress-concentration  fault  models  considered  in 


SPECTRAL 


Figure  14. 


Peak  slip  velocity  (low-passed,  5  Hz  cutoff)  and 
rupture  velocity  along  the  x  axis,  for  multiple 
stress-concentration  model.  Dashed  curves  for 
rupture  velocity  indicate  regions  which  ruptured 
out  of  sequence.  For  example,  rupture  occurred 
at  x  ■  5.6  km  while  the  region  between  5.0  and 
5.6  km  was  still  intact. 
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2.4  A  SIMULATION  STUDY  OF  THE  DETECTABILITY  OF  A  5.3  KT 
DECOUPLED  EXPLOSION  AT  REGIONAL  DISTANCES  IN  THE 
EASTERN  UNITED  STATES  BY  J.  R.  MURPHY  AND  T.  J.  BENNETT 

A  detailed  description  of  this  research  has  been  pre¬ 
viously  provided  in  the  Semiannual  Technical  Report  by  Bache, 
et  al.  (1980a) .  For  completeness,  we  include  below  an  execu¬ 
tive  summary  of  objectives  and  accomplishments. 


2.4.1  Objective 

The  objective  of  this  task  has  been  to  conduct  a  pre¬ 
liminary  theoretical  evaluation  of  CTBT  evasion  capabilities 
associated  with  the  cavity  decoupling  evasion  scenario.  In 
particular,  the  investigation  has  focused  on  extending  the 
existing  regional  seismic  data  base  for  explosions  in  the 
eastern  United  States  by  theoretically  scaling  observed  SALMON 
data  to  simulate  the  seismograms  to  be  expected  from  a  hypo¬ 
thetical,  fully  decoupled  5.3  KT  explosion  at  the  SALMON  shot- 
point. 

2.4.2  Accomp li shmen  t s 

Four  eastern  United  States  stations  were  selected  for 
analysis:  Eutaw,  Alabama  (EUAL,  A  =  246  km),  Cumberland 
Plateau  Observatory,  Tennessee  (CPO,  A  =  623  km) ,  Beckley, 

West  Virginia  (BLWV,  A  *  1065  km)  and  Wykoff,  Minnesota  (WFMN , 

A  »  1437  km) .  The  locations  of  these  four  stations  are  indi¬ 
cated  by  squares  on  Figure  17  which  provides  a  map  view  of  the 
North  American  Stations  which  recorded  the  SALMON  event  (Jordan, 
et  al. ,  1966).  The  observed  vertical  component  SALMON  data 
recorded  at  these  four  stations  have  been  theoretically  scaled 
(Murphy,  1977)  to  the  ground  motion  to  be  expected  from  a 
fully  decoupled  5.3  KT  explosion  at  the  SALMON  shotpoint.  The 
selected  theoretical  decoupling  factor  (with  respect  to  SALMON) , 
corresponding  to  a  simple  step  in  pressure  in  a  39  m  radius 


29 


SYSTEMS.  SCIENCE  AND  SOFTWARE 


cavity,  is  shown  in  Figure  18.  Observed  SALMON  vertical  com¬ 
ponent  seismograms  recorded  at  stations  EUAL  and  WFMN  are 
shown  at  the  top  of  Figures  19  and  20  respectively,  where 
they  are  compared  with  the  corresponding  simulated  decoupled 
seismograms  (bottom) .  The  simulated  decoupled  seismograms 
were  obtained  by  convolving  the  observed  SALMON  seismograms 
with  an  operator  defined  by  the  ratio  of  the  decoupled  to 
SALMON  theoretical  source  functions.  The  traces  shown  to 
the  left  of  the  simulated  decoupled  seismograms  correspond 
to  samples  of  the  noise  recorded  prior  to  SALMON  at  each 
station,  amplified  so  that  they  can  be  compared  with  the 
estimated  decoupled  signal  levels.  This  provides  a  basis 
for  assessing  the  detectability  of  the  hypothesized  decoupled 
explosion  at  these  stations  under  the  local  noise  conditions 
prevailing  at  the  time  of  the  SALMON  experiment.  A  review  of 
this  simulated  data,  as  well  as  that  from  the  other  two 
selected  stations,  indicates  that,  relative  to  the  noise  back¬ 
ground  prevailing  during  SALMON,  a  5.3  decoupled  explosion 
would  probably  have  been  detectable  on  the  short-period  re¬ 
cordings  from  the  two  nearer  stations  (EUAL  and  CPO)  but  not 
from  those  expected  at  the  more  distant  stations  (BLWV  and 
WFMN) .  With  respect  to  the  lower  noise  background  expected 
at  a  modern,  quiet  site,  the  Lg  phase  from  such  a  decoupled 
explosion  would  probably  be  detectable  at  all  four  stations, 
but  initial  P  wave  detection  would  still  be  questionable  at 
the  two  distant  stations.  Furthermore,  as  might  be  expected 
on  the  basis  of  previous  noise  studies,  increasing  the  effec¬ 
tive  bandwidth  of  the  short-period  data  to  encompass  lower 
frequencies  does  not  improve  this  detectability  due  to  the 
fact  that  the  noise  level  increases  with  decreasing  frequency 
while  the  decoupling  factor  remains  essentially  constant.  By 
the  same  argument,  consideration  of  higher  frequency  data 
should  improve  detectability,  at  least  for  high  Q  paths.  How¬ 
ever,  it  has  not  been  possible  to  quantitatively  evaluate 
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FREQUENCY  (Hz) 

Theoretical  decoupling  factor  for  5.3  KT 
in  a  39  m  radius  cavity  in  SALMON  salt. 
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Figure  19.  Comparison  of  observed  SALMON  and  simulated  decoupled  vertical 
component  seismograms.  Station  EUAL. 
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vigure  20.  Comparison  of  observed  SALMON  and  simulated  decoupled  vertical 
component  seismograms*  Station  WFMN. 


this  potential  improvement  using  the  available  data  set  due 
to  the  relatively  low  sampling  rate  employed  in  the  digitiza¬ 
tion. 
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III.  REGIONAL  SEISMOLOGY 


3.1  INTRODUCTION 

The  work  reported  in  this  section  pertains  to  Task  4.2. 
It  must  be  noted  that  several  of  the  numerical  calculations 
performed  under  Task  4.1.1  for  buried  explosions  have  been 
used  extensively  for  the  purposes  of  the  present  task  and 
have,  therefore,  not  been  described  in  Section  2.2  of  this 
report. 

The  results  pertinent  to  the  present  section  have  been 
the  object  of  a  detailed  S3  report  by  Bache,  et  al.  (1980b). 
They  have  also  been  summarized  in  the  Semiannual  Technical 
Report  by  Bache,  et  al.  (1980a) .  For  this  reason,  we  only 
include  for  completeness  an  executive  summary  in  the  following 
subsection. 


3.2  CONTRIBUTION  OF  TWO-DIMENSIONAL  SOURCE  EFFECTS  TO  THE 
FAR-FIELD  SEISMIC  SIGNATURES  OF  UNDERGROUND  NUCLEAR 
EXPLOSIONS  BY  T.  C.  BACHE,  T.  G.  BARKER,  N.  RIMER  AND 
J.  T.  CHERRY 


3.2.1  Objective 


As  the  contract  evolved,  the  objective  of  Task  4.2  be¬ 
came  to  study  the  far-field  seismic  waves  from  a  series  of 
two-dimensional  explosion  calculations  done  under  Task  4.1.1. 
The  purpose  was  to  determine  the  importance  of  nonlinear 
interaction  with  the  free  surface  (including  spallation)  and 
the  depth-dependence  of  overburden  pressure  and  material 
properties.  While  the  potential  importance  of  these  two- 
dimensional  effects  has  long  been  acknowledged,  previous 
attempts  to  model  the  seismic  signatures  of  underground 
explosions  have  almost  all  employed  spherically  symmetric 
point  source  representations.  A  quantitative  understanding 
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of  the  potential  influence  of  higher  order  source  effects  is 
important  for  delineating  the  validity  of  this  convenient 
approximation . 

3.2.2  Accomplishments 

The  results  summarized  here  are  extracted  from  a  1980 
S3  report  by  Bache,  Barker,  Rimer  and  Cherry.  The  report  pre¬ 
sents  a  detailed  analysis  of  the  seismic  waves  from  eleven  two- 
dimensional  finite  difference  calculations  of  underground 
nuclear  explosions  in  granite.  Seven  were  150  KT  explosions 
at  depths  from  159  to  1000  meters  done  by  J.  Trulio  and  N. 

Perl  of  Applied  Theory,  Inc.  (ATI). 

The  other  four  calculations  were  done  at  S3  with  one 
intended  to  model  the  PILEDRIVER  event.  The  constitutive 
model  and  source  geology  (three  layers)  were  chosen  for  this 
particular  event.  As  described  by  Rimer,  et  al.  (1979) ,  the 
computed  and  observed  ground  motions  were  compared  at  some 
twenty-five  near-field  gauge  locations,  where  the  agreement 
was  quite  good,  except  that  the  calculation  overpredicted 
the  amount  of  cracking  (spallatic.  This  turns  out  to  be 
important  when  comparing  synthetic  and  observed  far-field 
seismograms.  The  other  three  S3  calculations  were  the  same, 
except  that  the  depth  and  yield  were  varied.  The  four  calcu¬ 
lations  were  as  follows: 


DEPTH  YIELD 

(Meters)  (KT) 


463 

60  (PILEDRIVER) 

1000 

150 

1000 

20 

400 

20 
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The  main  results  of  the  study  may  be  summarized  by 
plots  of  Mg  and  m^  versus  source  depth  and  yield.  The  ob¬ 
served  values  for  the  United  States  granite  explosions 
PILEDRIVER,  SHOAL  and  HARDHAT  provide  a  standard  for  compari¬ 
son. 

In  Figure  21,  the  Mg  for  all  eleven  calculations  are 
plotted  versus  source  depth.  For  this  plot  all  M  ,  including 

9 

the  observed  value  for  PILEDRIVER,  have  been  scaled  to  150  KT 
by  adding  log  (150/W) . 

In  Figure  22,  the  Mg  values  are  plotted  versus  yield. 
Also  shown  are  the  Mg  from  seismograms  computed  with  reduced 
displacement  potential  (RDP)  sources.  The  RDP  calculations 
were  done  at  S3  for  two  depths  {460  and  1000  meters)  with  the 
same  constitutive  models  used  for  the  two-dimensional  calcula¬ 
tions.  Thus,  comparing  seismograms  from  the  RDP  and  the  two- 
dimensional  calculations  directly  displays  the  influence  of 
two-dimensional  effects. 

The  mb  data  from  the  ATI  and  S3  source  calculations  are 
plotted  versus  source  depth  in  Figure  23  and  versus  yield  in 
Figures  24  and  25.  Again,  the  values  for  comparable  RDP 
sources  are  shown  to  directly  display  the  two-dimensional 
effects  predicted  by  the  S3  calculations. 

The  most  important  result  of  this  study  is  that  the  S3 
calculations  show  Mg  to  be  a  strong  function  of  depth.  The 
shallow  S3  calculations  have  surface  wave  amplitudes  that  are 
a  factor  of  two  or  three  larger  than  those  from  a  comparable 
one-dimensional  source  calculation.  This  effect  is  probably 
exaggerated  because  the  free-surface  interaction  effects  are 
too  large  in  the  S3  calculations,  but  would  remain  important 
even  with  increased  strength  in  the  near-surface  material. 

On  the  other  hand,  the  ATI  calculations  show  no  strong  depen¬ 
dence  of  Mg  on  depth,  even  though  the  shallow  ATI  sources 
cratered.  We  do  not  know  why  the  two  sets  of  calculations  give 
such  different  results. 
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M„  VERSUS  LOG  YIELD 
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Figure  22.  The  Ms  values  for  the  S3  RDP  and  two-dimensional 
source  calculations  are  plotted  versus  explosion 
yield.  The  PILEDRIVER  observed  value  and  the 
range  of  Mg  for  the  seven  150  KT  ATI  sources  are 
also  shown.  The  source  depth  of  the  S3  two- 
dimensional  calculations  is  noted  with  the  Ms. 
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t  Figure  23.  The  ik  values  for  all  the  two-dimensional  calculations 

are  plotted  versus  source  depth.  The  S3  calculations 
at  the  same  yield  are  connected  with  a  dashed  line. 

The  observed  values  for  HARDBAT,  SHOAL  and  PILEDRIVER 
are  also  plotted  with  mg  *  -  0.25.  All  magnitudes 

have  been  scaled  to  150  KT. 
t 
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Figure  24.  The  b  phase  magnitude,  ra^,  is  plotted  versus  yield. 

The  values  for  the  ATI  calculations  fall  within  the 
indicated  range.  The  observed  values  are  the  pub-  | 

lished  m;:  minus  0.25  for  HARDHAT,  SHOAL  and  * 

PILEDRIVeR. 
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Figure  25.  The  c  phase  magnitude,  ra£  is  plotted  versus  yield. 

Published  (network  average)  magnitude  values  for 
HARDHAT,  SHOAL  and  PILEDRIVER  are  also  plotted. 
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Comparison  of  waveforms,  as  well  as  the  magnitude  data 
in  the  plots,  leads  to  the  following  conclusions: 


>lied  Theoi 


Inc.  Calculations 


•  Neither  mu  nor  Ms  are  strongly  dependent  on 
depth.  The  most  effect  was  on  Mg  at  shallow 
depths . 

•  Compared  to  the  semiempirical  R DP  source  of 
Mueller  and  Murphy  (1971) ,  depth  dependence 
in  both  amplitude  and  corner  frequency  is 
less  for  the  ATI  sources. 

•  The  pP  phase  appears  to  be  smaller  than 
expected  from  elastic  theory. 


Systems,  Science  and  Software  Calculations 

•  Two-dimensional  effects  are  not  very  important 
for  the  two  deep  explosions  (20  KT  and  150  KT 
at  lOOf'  meters)  .  Both  mb  and  Ms  are  little 
different  from  the  values  estimated  from  an 
RDP  source  computed  with  the  same  constitutive 
model  at  the  same  depth.*  It  is  gratifying 
that  these  two  very  different  and  complex  pro¬ 
cedures  arrive  at  the  same  results.  This  is 
true  even  though  considerable  cracking  and 
spallation  occur  in  the  two-dimensional  cal¬ 
culations  . 

•  For  body  waves  the  first  arriving  P  wave  is 
essentially  the  same  for  one-  and  two- 
dimensional  sources.* 

•The  two-dimensional  effects  enhance  the  sur¬ 
face  wave  amplitudes  for  the  shallow  events 
(60  KT  at  463  meters  and  20  KT  at  400  meters) 
by  a  factor  of  two  or  three.  We  must  quality 
this  by  pointing  out  that  these  shallow  cal¬ 
culations  have  very  strong  surface  interaction 
effects.  Comparison  with  PILEDRIVER  data  in¬ 
dicates  that  the  free  surface  interaction  is 


These  results  are  essentially  tests  of  the  entire 
computational  procedure. 
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overpredicted,  at  least  for  that  event.  Less 
free  surface  interaction  would  presumably  give 
less  enhancement  of  the  surface  wave  amplitudes. 

•  For  the  shallow  events,  the  m£  is  different 
than  predicted  with  an  RDP  source,  though  by 
less  than  1.2  units. 

•  Analysis  of  the  spectra  show  that  in  no  case 
is  pP  a  spectral  shadow  of  P,  as  it  is  for  an 
RDP  source  and  elastic  propagation. 

•  Phases  that  seem  to  be  associated  with  spall 
closure  can  be  seen  on  the  shallow  source  body 
wave  records.  However,  they  are  not  easily 
associated  with  identifiable  crack  closure 
patterns  in  the  source  calculation. 

•  Comparison  of  observed  and  calculated  body  and 
surface  waves  for  PILEDRIVER  leads  to  the  con¬ 
clusion  that  the  two  are  in  rather  good  agree¬ 
ment.  This  conclusion  is  subject  to  the  quali¬ 
fications  one  often  faces  in  this  kind  of  com¬ 
parison.  For  the  surface  waves,  it  is  the  need 
to  account  for  the  nonaxisywmetric  component, 
usually  attributed  to  tectonic  stress  release. 
Adding  a  recent  estimate  for  this  component  by 
Rivers  and  von  Seggern  (1979)  to  our  solution, 
we  get  good  agreement  with  the  data.  However, 
if  the  strike-slip  solution  of  Toksoz  and  Kehrer 
(1972)  is  correct,  the  synthetic  surface  waves 
are  about  a  factor  of  six  too  small. 

•  For  body  waves  the  comparison  is  complicated  by 
the  apparent  presence  of  strong  azimuthal  ef¬ 
fects  in  the  radiated  short  period  energy 
(Hadley  and  Hart,  1979) .  Our  conclusion  is 
that  the  computed  PILEDRIVER  source  (in  one-  or 
two-dimensions)  has  about  the  right  direct  P 
amplitude.  However,  the  two-dimensional  source 
calculation  appears  to  include  too  much  nonlinear 
interaction  with  the  free  surface.  The  later 
portion  of  the  P  waveform  does  not  match  the  data, 
apparently  because  pP  is  too  greatly  suppressed 
and  because  the  seismic  energy  from  spall  closure 
is  too  large  or  is  timed  incorrectly.  This 
"overprediction"  of  surface  interaction  effects 

is  expected  since  comparison  of  theoretical  and 
observed  near-field  motions  and  plots  of  the 
cracking  near  the  source  indicate  that  there 
was  too  much  spallation  in  the  calculation. 
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•The  constitutive  models  used  by  S3  (in  spherically 
symmetric  source  calculations)  lead  to  RDP  source 
functions  that  are  strongly  peaked,  with  the  value 
near  1  Hz  a  factor  of  five  or  more  larger  than  the 
value  at  long  periods.  The  peaking  is  due  to  the 
incorporation  of  an  effective  stress  law  and  the 
choice  of  unconfined  compressive  strength  (0.75  kb), 
based  on  laboratory  data  for  fractured  granite  and 
results  of  comparison  with  near-field  ground  motion 
observations. 

•  While  we  do  not  have  RDP  source  functions  for  the 
ATI  granite,  comparison  of  Ms  and  mt>  for  the  ATI 
two-dimensional  calculations  indicates  that  the 
RDP  peaking  is  probably  less  than  a  factor  of  two. 
Due  primarily  to  this  difference,  the  for  the 
S3  calculations  is  about  0.5  units  higher  than 
that  for  ATI  calculations  of  the  same  yield.  The 
ATI  Mg  values  fall  between  those  for  the  shallow 
and  deep  S 3  calculations . 


46 


SYSTEMS.  SCIENCE  AN O  SOETWAHE 


IV.  AUTOMATIC  SIGNAL  PROCESSING 


4.1  INTRODUCTION 

This  section  reviews  accomplishments  during  the  con¬ 
tract  period  on: 

TasJc  4.3.1  -  Development  of  the  MARS  Detector, 

Task  4.3.3  -  Review  of  Seismic  Event  Detectors. 

Although  we  had  hoped  initially  to  run  the  MARS 
(Multiple  Arrival  Recognition  System)  detector  at  VELA 
Seismological  Center  (VSC)  in  an  on-line  experiment  (Task 
4.3.2),  it  was  not  possible  to  perform  this  task.  The 
principal  reason  for  this  was  the  decision  to  commit  the 
resources  allocated  to  Task  4.3.2  to  the  MARS  detector 
task  (4.3.1).  The  results  which  were  obtained  from  this 
added  effort  were: 

1.  Increase  in  efficiency  obtained  by  develop¬ 
ment  of  heterodyne  technique. 

2.  Extensive  analysis  of  signal-to-noise  ratios 
using  the  modern  theory  of  statistical  sig¬ 
nal  analysis. 

Results  on  Task  4.3.1  are  described  in  the  next  sub¬ 
section,  where  the  MARS  detection  algorithm  is  described, 
together  with  an  evaluation  of  its  performance  in  an  off¬ 
line  mode  on  a  suite  of  test  tapes  supplied  by  VSC,  contain¬ 
ing  both  noise  and  signal  recordings.  A  more  complete  descrip¬ 
tion  may  be  found  in  Farrell,  et  al . ,  1980.  Results  on  Task 
4.3.3  have  been  detailed  in  a  separate  report  to  VSC  (Berger, 
1980)  and  we  give  below  an  executive  summary  of  the  most 
salient  points  and  recommendations. 


47 


srsrtMs.  science  ano  softwawc 


4.2  MARS  SEISMIC  DETECTION,  BY  W.  E.  FARRELL,  R.  C.  GOFF 

AND  J.  WANG 

4.2.1  Objective 

The  objective  of  this  task  was  to  optimize  the  MARS- 
based  signal  detection  algorithm  and  evaluate  it  in  an  off¬ 
line  mode  on  data  supplied  by  VSC. 

4.2.2  Accomplishments 

The  principal  accomplishment  has  been  the  development 
of  an  efficient  code  which  has  been  tested  on  four  hour  noise 
samples  from  NORSAR  and  Pinedale,  and  on  the  24-hour  NORSAR 
and  Pinedale  test  tapes  supplied  by  VSC.  Results  from  both 
test  tapes,  which  contained  over  30  signals  at  each  of  four 
different  amplitude  levels,  were  comparable  to  the  VSC  bench¬ 
mark  run;  that  is,  at  similar  false  alarm  rates,  about  the 
same  pattern  of  events  was  detected.  On  the  UNIVAC  1100/81 
the  algorithm  runs  at  40  times  real  time  on  the  10  sample/ 
second  NORSAR  data.  As  with  the  VSC  detector,  most  of  the 
time  (80  percent)  is  consumed  by  Fourier  transforms. 

As  an  example  of  the  current  state  of  the  MARS  detector. 
Figure  26  shows  the  first  ten  NORSAR  signals  extracted  from  the 
24-hour  test  tape.  Each  signal  was  buried  in  noise  at  four 
distinct  amplitude  levels,  which  Blandford  refers  to  as  "raw 
signal-to-noise  ratio  1/2,  1/4,  1/8,  1/16.”  We  refer  to  these 
four  scaled  copies  as  signal  classes  A,  B,  C  and  D.  Thus,  the 
plot  shows  ten  A  class  events,  those  with  raw  SNR  of  1/2. 
Furthermore,  the  seismograms  shown  in  this  figure  are  not  the 
actual  data  present  on  the  tape,  but  they  have  been  high-pass 
filtered  through  a  filter  with  a  1  Hz  corner  frequency.  This 
is  to  attenuate  the  long  period  energy  which  contains  mostly 
noise  power  and  very  little  signal  power.  The  events  are 
much  less  obvious  in  the  broadband  data  (see  Figure  29) .  For 
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each  seismogram  we  have  shown  the  times  at  which  "events" 
were  detected  by  the  MARS  processor.  Events  labeled  A  were 
detections  found  in  the  A  class  seismograms,  just  the  records 
shown  here.  Events  labeled  B  were  found  on  the  B  class  seismo¬ 
grams,  and  so  on. 

The  most  important  point  to  be  made  about  these  results 
is  that  the  principal  objective  in  the  work  to  date  has  been 
to  optimize  the  speed  of  the  code  and  to  calibrate  the  false 
alarm  behavior.  Thus,  we  have  not  "educated"  it  to  look  for 
specific  signal  types,  nor  to  pick  "first  motion."  Thus,  the 
narrow  band,  long  period  motion  in  Event  2,  is  ignored  in 
favor  of  later,  broadband  energy.  Seismogram  3  is  easily 
found  at  all  four  signal  levels,  but  the  peak  in  the  broad¬ 
band  power  falls  well  back  in  the  Lg  wave  train. 

It  should  be  emphasized  uhat  no  preprocessing  or  signal 
conditioning  (e.g. ,  optimum  band-pass  filtering)  precedes  the 
MARS  analysis.  Instead,  the  raw  data  is  at  the  very  first 
stage  decomposed  into  a  multiplicity  of  band  pass  channels. 

The  signal  amplitude  in  each  channel  changes  dynamically  as 
the  background  noise  power  rises  and  falls,  as  the  shape  of 
the  noise  spectrum  changes  with  time,  and  as  true  seismic 
events  come  and  go.  By  keeping  a  running  history  of  noise 
amplitude  in  each  channel,  the  MARS  process  then  flags  signifi¬ 
cant  short  term  excursions,  and  finally  searches  for  a  pattern 
across  the  many  frequency  bands. 

Among  the  subsidiary  accomplishments  in  this  automatic 
signal  processing  task,  the  most  important  are  the  following: 

•  A  tenfold  or  more  acceleration  in  processing 
speed  was  obtained  by  incorporating  a  hetero¬ 
dyne  and  decimate  operation  in  the  construction 
of  the  narrow  band  envelope  function.  Imple¬ 
mentation  of  this  feature  was  based  on  the 
observation  that  the  spectrum  of  a  narrow  band 
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signal  is  band  limited,  and  practically 
vanishes  outside  a  narrow  region  about  the 
center  frequency.  The  heterodyne  operation 
is  a  signal  processing  device  which  shifts 
the  center  frequency  down  to  zero  frequency. 

But  now  the  signal  spectrum  drops  to  a  negli¬ 
gible  value  at  frequencies  much  less  them 
the  Nyquist  frequency.  Then  the  inverse  trans¬ 
form  back  to  the  time  domain  can  be  accurately 
found  with  a  new  and  much  smaller  Nyquist  fre¬ 
quency  . 

•  An  extensive  study  was  made  of  the  statistical 
properties  of  narrow  band  seismic  noise.  This 
study  was  performed  on  the  continuous  envelope 
functions,  rather  than  the  envelope  maxima 
catalog  actually  used  by  MARS,  in  order  to  be 
able  to  relate  the  experimental  results  with 
classical  statistical  theory.  It  was  found, 
for  example,  that  the  envelope  probability 
density  function  closely  follows  the  Rayleigh 
law,  which  is  expected  if  narrow  band  envelopes 
are  extracted  from  Gaussian  noise  whose  spectrum 
is  not  too  ragged.  Further,  the  joint  statis¬ 
tics  of  multiple  envelope  functions  were  studied 
to  estimate  the  cross-correlation  between  two 
envelope  functions  generated  by  different  filters. 
This  work  was  necessary  to  quantify  the  possi¬ 
bility  of  cross-correlation  in  envelope  maxima, 
and  to  give  objective  guidance  both  to  the  con¬ 
struction  of  the  comb  of  filters,  and  to  the 
significances  of  envelope  alignments.  Preliminary 
work  was  undertaken  in  multiple  envelope  correla¬ 
tion,  using  the  Bernoulli  trial  model,  but  we  have 
not  yet  obtained  definitive  results. 
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•  Working  still  with  noise  records,  we  know  how  to 
tune  the  MARS  false  alarm  rate  by  adjusting  the 
three  key  parameters:  amplitude  (normalized 
with  respect  to  the  noise  in  each  band) ,  band¬ 
width  and  dispersion.  Figure  27  shows  a  cube 

in  this  parameter  space,  each  cell  of  which 
contains  a  number  related  to  the  False  Alarm 
Rate.  The  bigger  this  number ,  the  poorer  the 
detector.  In  our  terminology,  the  bandwidth 
parameter  is  called  k;  it  gives  the  number  of 
envelope  maxima  (out  of  twenty  possible)  re¬ 
quired  before  an  event  is  called.  The  amplitude 
parameter,  also  referred  to  as  normalized  ampli¬ 
tude  cutoff  in  the  third  figure,  expresses  how 
large  extrema  must  be  before  an  envelope  peak 
is  to  be  included  in  the  search  for  events. 

We  refer  to  this  as  a  normalized  amplitude  be¬ 
cause  each  envelope  function  in  the  time  window 
currently  being  processed  is  scaled  by  the  mean 
envelope  amplitude  of  the  previous  N  windows. 

There  is,  then,  a  floor  on  the  minimum  instan¬ 
taneous  signal  spectrum  required  before  an 
event  is  called.  The  floor  is  just  some  multiple 
of  the  preceding  noise  spectrum.  Finally,  the 
dispersion  parameter  specifies  how  broad  a  time 
swath  one  permits  the  k  peaks  to  fall  within. 

•  Figure  28  shows  a  specific  result  selected  from 

our  noise  study.  It  is  a  slice  along  the  k  «  6  ; 

plane  for  a  two-hour  section  of  NORSAR  data.  i 

The  figure  shows  that  the  lower  left  hand  corner  ; 

is  good,  but  that  the  upper  right  hand  corner  ■ 

is  bad.  But  we  do  not  want  to  be  too  far  down 
in  the  lower  left  hand  corner.  That  gives  great 
noise  rejection,  but  concomitantly  requires  an 
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Figure  27.  Schematic  representation  of  the  three-dimensional  j 

parameter  space  spanned  by  the  key  MARS  detection  i 

parameters.  Each  cell  within  this  cube  contains  j 

a  number  giving  the  false  alarm  rate  when  the  ! 

algorithm  operates  on  pure  noise.  j 
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extremely  large  and  impulsive  arrival  to  trigger. 
Thus,  this  part  of  the  parameter  space  would  show 
a  low  probability  of  detection  for  easily  seen 
signals  To  capture  weak,  moderately  dispersed 
signals,  one  wants  to  site  the  detector  param¬ 
eters  somewhere  in  the  hatched  region. 

e  New  graphical  display  procedures  have  been  imple¬ 
mented  for  showing  the  time-frequency  structure 
of  se i smog r ams .  These  have  been  found  to  be  ex¬ 
tremely  powerful  tools  in  our  research  studies, 
and  should  be  equally  useful  at  seismic  analyst 
stations  for  operational  detection  and  discrimina¬ 
tion.  Figure  29  displays  the  prebandpass  filtered 
version  of  Seismogram  5  from  Figure  26,  along 
with  its  associated  time- frequency  decomposition. 
This  is  a  true  amplitude  plot  at  each  frequency, 
which  clearly  indicates  the  dominance  of  the 
long  period  noise,  as  does  the  comparison  between 
this  version  of  the  seismogram  and  the  bandpassed 
version  in  Figure  26.  It  can  be  seen,  for  ex¬ 
ample,  that  the  best  signal  amplitude  occurs 
over  the  2.0  to  3.5  Hz  frequency  band.  There 
is  also  significant  signal  power  out  to  5.0  Hz, 
but  the  noise  there  is  so  low,  it  is  not  possible 
to  get  an  accurate  estimate  of  the  signal-to- 
noise  ratio. 

•  Another  type  of  display  is  shown  in  Figure  30. 

This  display  is  similar  to  the  sonogram  display, 
long  used  in  acoustic  signal  processing,  but 
these  figures  were  generated  by  digital  rather 
than  analog  signal  analysis  methods.  In  this 
display,  high  power  points  in  the  time-frequency 
plane  are  shown  as  black  zones,  small  amplitude 
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Figure  29.  Seismic  Event  S  (see  Figure  26)  and  its  associated 
time- frequency  decomposition.  The  fin  at  16  sec 
(this  time  is  arbitrary)  and  2.5  Hz  is  the  same 
event  which  appears  at  56  seconds,  2.5  Hz  in  the 
middle  panel  of  Figure  30.  In  this  figure  the  true 
amplitude  in  each  frequency  band  is  preserved. 

This  makes  the  predominance  of  the  long  period 
noise  obvious.  In  Figure  30,  the  data  has  been 
"whitened"  by  scaling  according  to  the  largest 
amplitude  in  each  band.  This  emphasizes  frequency 
bands  with  good  signal- to-noise  ratio. 
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Figure  30.  "Sonograms"  for  seismic  Events  3,  5  and  9  which 
are  also  displayed  in  Figure  26. 
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regions  shade  down  through  various  gray  levels 
towards  pure  white.  Unlike  the  previous  dis¬ 
play,  these  data  have  been  scaled  according  to 
the  "peak"  power  in  each  frequency  band.  This 
tends  to  "whiten"  any  signals  present,  and  make 
more  visible  these  frequency  bands  with  good 
signal- to-noise  ratio. 

•  The  extra  frequency  dimension  in  Figures  29  and 
30  tells  much  more  about  the  nature  of  the 
seismic  events  than  does  the  seismogram  itself. 
Event  9,  for  example,  is  barely  visible  in 
Figure  26,  but  on  Figure  30,  it  shows  up  as 
the  dark  band  between  2.0  and  3.0  Hz,  slightly 
later  than  nine  minutes.  One  does  see,  of  course, 
bands  of  similar  length  elsewhere  in  the  figure, 
but  they  are  not  surrounded  by  the  long  gray  zone, 
indicating  they  are  most  likely  noise  signals. 


4.3  DETECTION  STATE-OF-THE-ART  BY  J.  BERGER 
4.3.1  Objective 

The  objective  of  this  review  study  was  to  identify 
all  pertinent  automatic  signal  detection  algorithms  by  equa 
tion,  mathematical  derivation  and  usage  history,  including 
algorithms  currently  under  development  and  scheduled  to  be 
validated  by  October  1981. 


4.3.2  Accomplishments 

A  draft  report  entitled  "Seismic  Detectors"  The 
State-of-the  Art"  has  been  prepared  and  submitted  to  VSC. 
In  this  report,  the  following  detection  algorithms  were 
reviewed: 
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1. 


IBM 


2.  Z 

3.  Deflection 

4.  Analytic  Envelope 

5.  Allen 

6.  Stewart 

7.  Walsh 

8.  MARS 

General  conclusions  reached  were: 

1.  All  detectors  reviewed  in  this  report  are  based  on 
a  comparison  between  some  form  of  variance  estimate  of  the 
signal  calculated  over  a  time  period  approximately  equal  to 
the  duration  of  an  expected  seismic  event,  and  the  normal 

(or  long  term)  variance. 

Possible  exceptions  to  this  general  form  are  the 
phase  detector  described  by  Unger  (1978)  and  the  similar  part 
of  MARS .  The  former  did  not  work  because  the  short  period 
phase  fluctuations  were  so  rapid  that  they  could  not  reliably 
be  predicted;  the  latter  has  not  been  implemented. 

2.  Ad  hoc  detection  algorithms  such  as  Allen's  and 
Stewart's  (1977)  offer  a  hope  of  some  improvement  in  the 
detection  process.  These  methods,  after  having  made  the 
basic  detection  statistic  test  and  having  made  a  tentative 
identification,  are  followed  by  a  further  process  that  is 
designed  to: 

a.  reduce  the  false  alarm  rate, 

b.  improve  the  timing  capability,  and 

c.  speed  the  recovery  of  the  detection 
algorithm  after  an  event  is  encountered. 
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How  much  these  algorithms  can  be  expected  to  offer 
significant  improvements  depends  upon  the  noise  environment 
in  which  they  are  operating.  If  the  SNR  is  such  that  a  high 
probability  of  detection  is  achieved  with  a  tolerable  FAR, 
then  little  improvement  can  be  expected.  If,  however,  the 
SNR  is  poor,  so  that  the  probability  of  detection  is  low  for 
a  given  FAR,  then  significant  improvement  may  be  achieved  by 
such  two-stage  algorithms. 

3.  Testing  of  the  detectors  described  in  this  report 
was  not  carried  out  in  a  uniform  manner  on  the  same  or  even 
similar  data  sets.  Thus,  no  definitive  statement  about  their 
relative  performance  may  be  made.  However,  it  appears  that 
no  one  detector  is  obviously  superior  to  all  others.  Some 
have  been  optimized  for  teleseismic  detection  using  a  moder¬ 
ately  large  computing  facility.  Some  have  been  designed  for 
operation  in  microprocessor  controlled  field  recording  units, 
typically  used  in  the  near  field  or  at  least  at  regional 
distances.  These  application  stress  ease  of  implementation 
in  microprocessors.  Some  algorithms  stress  timing  ability  as 
this  may  be  another  object  of  the  algorithm. 

4.  Little  theoretical  advance  has  been  achieved  since 
Frieberger's  1963  work  in  which  he  "solved"  the  problem  for 
detection  of  a  Gaussian  signal  in  Gaussian  noise.  The  influ¬ 
ence  of  this  work  has  rightly  guided  the  design  and  implementa¬ 
tion  of  many  of  the  detectors  used  today.  However,  it  must 
not  be  thought  that  the  "real  problem"  has  been  solved.  The 
"real  problem"  is  the  detection  of  certain  nonGaussian  signals 
in  the  presence  of  nonGaussian  noise.  Frieberger  discovered 
the  optimum  detector  for  an  approximate  model  of  the  actual 
situation,  but,  other  algorithms  may  work  significantly  better 
on  "real"  data.  Further,  considerations  such  as  timing  ability, 
recovery  of  the  algorithm  after  an  event,  and  immunity  to 
highly  nonGaussian  noise  such  as  line  spikes  and  data  dropouts 
may  dictate  very  different  approaches. 
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As  a  result  of  the  review,  the  following  recommendations 
were  made: 

•  Uniform  testing  of  all  viable  detectors  should 
be  conducted  with 

*  a.  a  realistic  synthetic  data  set; 

b.  all  detectors  coded  for  and  running 
on  the  same  or  similar  machines; 

c.  ROC  curves  produced  for  each  detector 
comparing  them  with  the  theoretical 

*  performance  of  the  detector  for  a  com¬ 
pletely  known  signal. 

•  Theoretical  and  experimental  research  in  phase 
sensitive  detectors  should  be  supported  as  this 

I  is  an  area  where  it  may  be  possible  to  exploit 

signal  information  not  used  in  energy  or  power 
detectors. 

•  Research  should  be  encouraged  to  test  if  matched 
filters  for  specific  types  of  signals  (station/ 

■  source  pairs)  improves  detector  performance. 

•Further  development  of  hybrid  detection  algorithms, 
that  combine  high  probability  of  detection  with  a 
high  FAR  and  then  are  followed  by  a  post  processing 

i  to  reduce  the  FAR  should  be  undertaken. 

•  New  algorithms,  designed  specifically  for  use 
with  the  new  generation  of  three  component  broad 
band  seismographs,  should  be  developed  and  tested. 

•  In  view  of  data  rates  from  these  new  data  sources, 
dedicated  microprocessors  to  run  the  detection 
algorithms  on  a  one  processor  per  channel  (or 
seismic  station)  basis  should  be  considered. 
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V.  EVALUATION  OF  SUITES  OF  EVENTS 


5.1  INTRODUCTION 

The  objective  of  this  task  (Task  4.4  of  our  Contract 
F08606-79-C-0008)  was  to  study  and  evaluate  the  feasibility 
of  automatic  determination  of  body  wave  magnitude  using  the 
MARS  processor  described  in  Section  4.2.  In  addition,  the 
work  included  developing  and  testing  of  an  analogous  auto¬ 
matic  surface  wave  magnitude  algorithm. 

This  work  is  summarized  in  an  S3  Semiannual  Technical 
Report  (Bache,  et  al. ,  1980a).  The  next  subsection  summarizes 
the  most  salient  results. 


5.2  AUTOMATED  MAGNITUDE  in.  AND  M  ,  BY  T.  C.  BACHE  AND 

B.  SHKOLLER 

In  FY  *79,  we  proposed  an  automated  magnitude  measure 
called  mfa  and  tested  it  by  processing  high  signal/noise  re¬ 
cordings  of  eleven  Pahute  Mesa  explosions  at  six  teleseismic 
stations  (Bache,  1979;  Bache,  Day  and  Savino,  1979) .  We  con¬ 
cluded  that  the  mb  was  at  least  as  good  a  magnitude  measure 
and  yield  indicator  as  the  most  carefully  determined  time 
domain  m^.  The  m^  algorithm  used  in  that  earlier  work  did 
not  account  for  the  presence  of  seismic  noise,  and  so  was 
viewed  as  primarily  a  demonstration  of  the  concept. 

As  part  of  the  VSC  conducted  discrimination  experiment, 
algorithms  were  constructed  to  compute  variable  frequency 
magnitudes,  m^tf),  at  selected  frequencies  (Masso,  et  al. , 
1979;  Savino,  et  al. ,  1980).  These  algorithms  do  incorporate 
corrections  for  the  presence  of  interfering  phases  and  seismic 
noise.  Otherwise  the  mb  used  by  Bache  (1979)  is  not  much 
different  from  mb  (1  Hz) . 
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In  the  current  contract,  we  have  implemented  a  new 
algorithm  for  (which  is  essentially  (1  Hz))  and  tested 
it  by  application  to  a  series  of  synthetic  seismograms  which 
are  constructed  by  embedding  a  signal  of  known  m^  in  seismic 
noise.  The  m^  has  also  been  computed  for  a  set  of  RKON  re¬ 
cordings  of  Eurasian  events  taken  from  the  AI  data  set.  A 
similar  algorithm  has  been  implemented  for  computing  an 
analogous  M  ,  and  some  testing  has  been  done. 

S 


5.2.1  Description  of  the  m^  Algorithm 


The  calculation  of  m^  includes  the  following  steps: 

1.  The  seismogram  is  processed  by  MASS,  which 
applies  a  suite  of  narrow  band  filters  and 
constructs  a  table  of  envelope  function 
peak  amplitudes  (A  )  and  group  arrival  times 

/  a.  \  ^ 


2.  The  MARS  detection  algorithm  is  applied  to 
identify  an  undispersed  P  wave  arrival. 

This  algorithm  works  with  a  "surface"  like 
that  is  shown  in  Figure  31,  where  the  Ac 
are  plotted  versus  tg  and  filter  center0 
frequency,  fg.  For  lach  fc  the  largest 

Ac  is  denoted  by  an  asterisk.  The  relative 
amplitudes  of  the  other  peaks  are  denoted 
by  integers,  with  the  asterisk  given  the 
value  10.  The  detection  algorithm  searches 
for  a  "ridge"  of  amplitudes  with  nearly  the 
same  tg  to  find  the  group  arrival  time,  €g, 
of  a  detected  P  wave.  The  "+"  and  "-" 
symbols  in  the  plot  denote  a  statistically 
determined  uncertainty  window  within  which 
any  Ac  could  be  associated  with  the  detected 
signal. 

3.  For  each  fc,  one  Ac  is  chosen  to  define  the 
"MARS  signal  spectrum. "  This  is  done  by  a 
scoring  algorithm  which,  in  its  current  form, 
works  fairly  well,  but  could  probably  be  im¬ 
proved  for  low  S/N  events. 
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Figure  31. 


10  sac 


Bottom  portion  is  a  seismogram  recorded  at  a  station 
in  Bluff,  Alaska,  from  an  event  in  the  Kurils.  The 
top  half  is  a  plot  of  filter  center  frequency  versus 
group  arrival  time  (sec) .  The  t|  with  arrow  attached 
shows  where  the  automated  detection  and  first  arrival 
algorithms  picked  the  signal  arrival  time. 
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4.  The  Ac  values  about  a  specific  fc  (e.g., 

1  Qz)  are  fit  to  give  a  "smoothed"  esti- 
ment  for  the  spectral  amplitude  at  that 
frequency.  This  amplitude  is  corrected 
for  the  presence  of  nearby  envelope  peaks 
and  a  statistical  estimate  of  the  uncer¬ 
tainty  is  made,  based  on  the  "MARS  spectrum" 
of  the  ambient  noise.  The  m^  is  computed 
from  log  (A~**;  fc)  +  where  Ag**  is  the 
noise  corrected  Ac  and  B  is  the  Usual  dis¬ 
tance  correction. 

5.2.2  m^  Algorithm  Testing 

In  Figure  32,  we  show  the  for  a  synthetic  seismogram 
which  has  been  added  to  a  noise  sample  from  RKON,  with  the 
ratio  of  peak  time  domain  amplitudes  ("Peak  S/N")  being  scaled 
to  100.  The  mb  is  3.589,  which  may  be  compared  to  the  time 
domain  of  3.78,  computed  from  log(A/T)  +  B,  with  T  =  0.7 
seconds.  Systematic  differences  between  nv  and  m^  are  ex¬ 
pected. 

The  same  synthetic  was  added  to  various  segments  of 
noise  recorded  at  RKON,  Some  examples  are  shown  in  Figure  33. 
Bache,  et  al.  (1980)  show  fifteen  such  examples.  The  error  in 
arrival  time  and  for  these  examples  are  tabulated  in  Table 
1.  Only  when  the  ratio  of  MARS  signal  spectrum  to  MARS  noise 
spectrum  gets  to  be  1.5  or  smaller  are  the  errors  very  large. 
That  is,  even  when  the  signals  can  barely  be  detected  by  eye, 
the  algorithm  is  able  to  give  a  good  estimate  for  the  mag¬ 
nitude. 

5.2.3  for  Actual  Data 

Bache,  et  al.  (1980)  present  the  m^  results  for  Eurasian 
explosions  and  earthquakes  recorded  at  RKON  and  LASA.  The  epi- 
central  distances  are  from  53  to  94  degrees  and  the  events  have 
PDE  magnitudes  from  4.8  to  6.0.  The  algorithm  appears  to 
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Figure  32. 


Calculation  of  for  a  seismogram  with  peak  synthetic/peak 
noise  *  100.  The  seismogram  is  shown  at  the  top  with  an 
asterisk  denoting  the  signal  arrival  time  tg.  The  mv  calcu¬ 
lation  is  shown  graphically,  with  the  heaviest  line  being  log 
(A£;  fc)  +  B.  The  line  of  medium  shade  is  the  analogous  mag¬ 
nitude  for  the  noise.  The  light  line  is  the  tg  for  each 
A|*(fc) ,  with  zero  on  the  tg  scale  (at  right)  Being  tg.  The 
mb  is  computed  at  1  Hz  and  Is  marked  with  an  X.  ^Printed 
above  the  graph  are  mb  (MB) ,  the  uncertainty  in  mb  imposed 
by  the  noise  (+  DMB),  the  noise  magnitude  at  1  Hz  (MBNF) , 
the  mean  signaT/noise  power  over  the  frequency  band  from 
FLEFT  to  FRIGHT,  and  the  B  correction  for  the  magnitude 
(BDEL) . 
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Figure  33.  Test  of  the  mv,  algorithm  with  noise  Section  62110.  The  ratio  of 
the  peak  amplitude  of  the  synthetic  to  the  peak  noise  amplitude 
is  denoted  by  Peak  S/N.  The  actual  arrival  time  of  the  signal 
is  the  time  of  the  asterisk  plus  t  . 


TABLE  1 


SUMMARY  OF  ^  TEST  CALCULATIONS 


Identifier 

Peak  Siqnal 
Peak  Noise 

S/N 
(1  Hz) 

“b 

”b 

Uncertainty 

"b 

Error 

Arrival 

Time 

Error 

Noise 

Section 

62110 

(Figure  32) 

100 

217 

3.589 

0.001 

- 

— 

A 

3 

6.5 

3.583 

0.018 

-0.006 

+0.02 

B 

1 

2.2 

3.593 

0.055 

0.004 

0.03 

C 

1 

2.6 

3.667 

0.046 

0.078 

-0.09 

D 

1 

2.2 

3.586 

0.056 

-0.003 

0.48 

E 

1 

2.2 

3.581 

0.057 

-0.008 

-0.01 

F 

0.5 

1.5 

3.734 

0.082 

0.145 

0.09 

G 

0.5 

1.2 

3.622 

0.110 

0 . 03  j 

-0.76 

H 

0.5 

1.5 

3.708 

0.088 

0.119 

0.14 

I 

0.5 

1.4 

3.689 

Missed  Signal 
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give  excellent  results,  even  when  the  signal  is  difficult  to 
detect  by  eye  so  that  time  domain  m^  measurements  would  prob¬ 
ably  not  be  attempted. 


5.2.4  Description  of  the  Mg  Algorithm 

The  Mg  algorithm  is  nearly  the  same  as  the  m^  algorithm, 
except  that  the  dispersion  character  of  the  surface  wave  sig¬ 
nal  of  interest  must  be  specified  explicitly.  The  t  -fc  plane 
(Figure  31)  is  first  mapped  into  a  t'-f  plane  where  a  signal 

y  C 

with  the  expected  dispersion  appears  to  be  nearly  undispersed. 
The  algorithm  then  follows  the  same  steps  as  for  m^. 

5.2.5  Mg  Algorithm  Testing 

In  Figure  34,  we  show  an  example  of  the  Mg  calculation. 

A  synthetic  seismogram  was  added  to  a  long  period  noise  segment 
(from  LASA) ,  with  the  peak  signal  and  noise  amplitudes  being 
equal.  Comparing  to  the  Mg  calculation  for  the  synthetic  it¬ 
self,  the  20  second  arrival  time  error  is  4  seconds  and  the 

A 

error  is  0.06.  This  and  similar  calculations  done  with  a 
s 

"perfect"  estimate  for  the  dispersion  of  the  expected  signal, 

A 

show  that  the  Mg  algorithm  works  very  well.  Testing  with  real 
data  where  the  prespecified  dispersion  is  only  approximately 
correct  has  not  yet  been  done.  However,  we  point  out  that  the 
Mg  calculation  only  requires  an  accurate  identification  of  the 
MARS  spectrum  in  a  narrow  period  band,  in  this  case  near  20 
seconds. 

5.2.6  Conclusions 

The  automated  magnitude  algorithms  appear  to  give  con¬ 
sistent  and  convenient  results,  though  much  more  extensive 
testing  with  real  data  is  required.  The  and  Mg  values  may 
be  viewed  as  a  by-product  of  signal  processing  with  the  MARS 
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program,  which  is  also  configured  for  signal  detection  and 
earthquake/explosion  discrimination.  Thus,  these  values  can 
easily  be  catalogued  when  processing  data  with  MARS  for  other 
purposes. 
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VI.  GROUND  MOTION  ANALYSIS 


6 . 1  INTRODUCTION 

Task  4.5  involved  a  compilation  and  assessment  of  near 
source  data  recorded  from  underground  explosions  with  the  pur¬ 
pose  of  determining  their  significance  for  the  definition  of 
nuclear  source  functions  for  different  emplacement  media.  More 
specifically,  previously  unpublished  data  from  the  Sandia 
Corporation  archives  are  being  reviewed  and  integrated  into 
the  previously  compiled  data  base. 

Work  on  this  task  was  the  topic  of  a  report  by  Murphy 
and  Bennett  (1980)  and  was  described  by  Bache,  et  al.  (1980a) . 

A  summary  is  provided  in  the  following  subsection. 


6.2  ANALYSIS  OF  FREE-FIELD  DATA  FROM  EXPLOSIONS  IN  WET 

TUFF  AND  ALLUVIUM  EMPLACEMENT  MEDIA  BY  J.  R.  MURPHY 

AND  T.  J.  BENNETT 

Table  2  lists  the  events  for  which  previously  unpub¬ 
lished  free-field  data  have  been  collected.  It  can  be  seen 
that  data  from  eleven  different  explosions  covering  a  yield 
range  from  about  1  to  1200  KT  are  being  investigated.  These 
events  were  detonated  in  wet  tuff /rhyolite  and  alluvium  emplace¬ 
ment  media  and  the  associated  depths  of  burial  vary  by  more  than 
a  factor  of  six.  These  data,  together  with  those  previously 
collected,  are  being  used  to  assess  the  dependence  of  seismic 
coupling  on  source  medium.  For  example,  Figure  35  shows  a 
comparison  of  the  reduced  displacement  potentials  (RDP's)  pre¬ 
dicted  using  an  analytical  model  based  on  Pahute  Mesa  explo¬ 
sions  in  wet  tuff /rhyolite  with  several  observed  RDP's  from 
Yucca  Flat  events  in  wet  tuff.  It  can  be  seen  that  the  pre¬ 
dicted  and  observed  agree  fairly  well,  which  suggests  that 
the  seismic  coupling  properties  of  Yucca  Flat  tuff  are  similar 
to  those  of  the  Pahute  Mesa  volcanics.  Howev-  . ,  a  very  different 
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TABLE  2 


NUCLEAR  EXPLOSIONS  FOR  WHICH  PREVIOUSLY  UNPUBLISHED 
FREE-FIELD  DATA  HAVF.  BEEN  COLLECTED 


EVENT 

MEDIUM 

YIELD 

DEPTH 

AGILE 

Wet  Tuff 

L-I 

732  m 

BOXCAR 

Wet  Rhyolite 

1200  KT 

1165  m 

CALABASH 

Wet  Tuff 

110  KT 

626  m 

CARPETBAG 

Wet  Tuff 

220  KT 

663  m 

COMMODORE 

Wet  Tuff 

250  KT 

749  m 

CYPRESS 

Wet  Tuff 

L 

411  m 

DIESEL  TRAIN 

Wet  Tuff 

L 

419  a 

FAULTLESS 

Wet  Tuff 

I 

975  m 

HUTCH 

Wet  Tuff 

I 

549  m 

LANPHER 

Wet  Tuff 

L-I 

714  m 

PLATTE 

Tuff 

1.85  KT 

191  m 
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PREDICTED  WET  T/R 


picture  emerges  when  this  same  procedure  is  applied  to  explo¬ 
sions  in  dry  Yucca  Flat  alluvium.  Figure  36  shows  a  compari¬ 
son  of  the  observed  RDP's  from  four  events  in  dry  alluvium 
with  the  corresponding  RDP's  predicted  for  the  same  yields 
and  depths  of  burial  in  wet  tuff/rhyolite  emplacement  media. 

It  can  be  seen  that  these  data  clearly  indicate  that  dry 
alluvium  is  a  weak  coupling  medium  with  respect  to  wet  tuff/ 
rhyolite,  in  agreement  with  the  observed  differences  in  the 
teleseismic  mfa  data  from  explosions  in  these  two  media.  Thus, 
the  available  free-field  data  are  essentially  consistent  with 
the  teleseismic  m^  data  and,  consequently,  should  be  of  use 
in  their  interpretation. 
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VII.  ACOUSTIC  WAVES  FROM  EXPLOSIONS  IN  THE  OCEAN 

7.1  INTRODUCTION 

The  study  o£  the  acoustic  signature  of  a  near  surface 
oceanic  explosion  gave  rise  to  a  report  by  Bache,  et  al.  (1980c) . 
An  executive  summary  of  the  results  is  included  below. 


7.2  THE  UNDERWATER  ACOUSTIC  SIGNATURE  OF  A  NUCLEAR 
EXPLOSION  AT  THE  OCEAN  SURFACE  BY  T.  C.  BACHE, 
T.  G.  BARKER,  M.  G.  BROWN,  K.  D.  PYATT  AND 
H.  J.  SWANGER 


7.2.1  Objective 


The  objective  of  this  task  is  to  estimate  the  gross 
spectral  character  and  duration  of  the  acoustic  wave  signa¬ 
ture  of  a  nuclear  explosion  near  the  ocean  surface.  The 
explosion  is  assumed  to  have  a  yield  of  1  KT  and  the  nominal 
range  is  6600  km. 


7.2.2  Accomplishments 


The  main  results  of  this  work  are  summarized  in  a 


report  by  Bache,  et  al.  (1980) .  The  technique  was  to  con¬ 
struct  theoretical  pressure-time  histories  by  convolving  am 
elastic  representation  for  the  source  with  a  Green’s  function 
representing  the  propagation.  The  source  may  include  two 
parts.  The  first  is  the  airblast  loading  applied  to  the  sur¬ 
face  and  the  second  is  the  direct  coupling  of  energy  into  the 
water.  When  the  explosion  is  above  a  certain  height,  esti¬ 
mated  to  be  2.3  meters  for  1  KT,  the  direct  coupling  contri¬ 
bution  may  be  ignored.  Thus,  the  calculations  include  the 
cases  when  the  explosion  is  on  or  above  the  surface,  though 
no  attempt  is  made  to  account  for  height-of-burst  effects  that 
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may  become  important  when  the  (1  KT)  explosion  is  more  them 
10  meters  above  the  surface. 

Both  components  of  the  source  are  expressed  in  terms 
of  pressure- time  histories,  P(r,t)  applied  to  the  ocean  sur¬ 
face.  The  total  force  applied  is 

roo 

F(id)  -  j  P(r,u)  JQ(ojpr)  rdr  .  (3) 

■'o 

The  spectral  amplitudes  of  the  total  force  for  both  components 
of  the  source  are  plotted  in  Figure  37  for  the  ray  parameter 
(p  *  0.667  sec/km)  appropriate  for  large  distances. 

The  pressure  at  the  receiver  is  obtained  by  convolving 
F(o>)  with  the  appropriate  Green's  function,  which  is  the  pres¬ 
sure  at  the  receiver  due  to  a  vertical  point  load  applied  at 
the  center  of  the  source  region.  The  Green's  functions  were 
computed  with  the  WKBJ  seismogram  method  (Chapman,  1978; 

Brown,  et  al. ,  1980) .  The  results  turned  out  to  be  very  little 
different  from  those  from  geometric  ray  theory.  Then  for  any 
particular  ray  arrival,  the  spectral  pressure  at  the  receiver 
is  approximately 

Pr  (oj  )  a  i  w  F(uj)  •  Q(ui)  ,  (4) 

where  Q(oj)  represents  the  frequency-dependent  filter  imposed 
by  the  travel  path.  However,  at  the  frequencies  of  interest 
(<  50  Hz),  absorption  in  the  ocean  is  very  small.  Therefore, 
we  assume  that  Q(u»)  -  1.  The  total  record  is  a  summation  of 
many  rays,  leading  to  a  total  spectrum  which  is  merely  a 
modulated  version  of  Equation  (4) . 

In  Figure  38,  we  plot  our  estimates  for  the  spectra  of 
the  pressure  pulse  at  a  hydrophone  in  the  deep  ocean  at  a 
large  distance  from  the  explosion.  The  direct  coupling 
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FREQUENCY  (101  Hz) 


Figure  37.  The  explosion  source  spectra. 
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PRESSURE  (ARBITRARY  UNITS) 


FREQUENCY  (Hz) 


Figure  38.  The  spectrum  of  the  pressure  time  history  of  a  single 
ray  arrival  at  a  SOFAR  axis  receiver  station.  The 
complete  signal  spectrum,  which  is  a  sum  of  many  single 
ray  spectra,  would  have  the  same  basic  shape.  The 
amplitude  units  are  arbitrary,  depending  on  the  amplitude 
of  the  Green's  function  for  the  travel  path. 
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contribution  is  small  and  the  spectrum  is  band-limited,  having 
a  peak  near  22  Hz  and  being  down  by  a  factor  of  3  from  the 
peak  at  about  4  and  50  Hz.  The  spectra  cube-root  scale  with 
yield.  That  is,  if  the  yield  were  2  KT,  the  spectrum  in 
Figure  38  shifts  to  frequencies  smaller  by  a  factor  of  1.26. 

To  estimate  the  signal  duration,  the  Green's  functions 
were  computed  with  WKBJ  seismogram  theory  for  a  variety  of 
assumptions  about  the  ocean  sound  velocity  profile.  At  the 
large  ranges  (>  3500  km)  of  interest,  the  main  contribution 
is  from  rays  that  turn  above  the  ocean  bottom.  For  a  reflec¬ 
tion  coefficient  of  as  large  as  0.98,  the  bottom  bounces  make 
almost  no  contribution  because  they  reflect  so  many  times. 

The  calculated  signal  duration  varies  widely,  depending 
upon  the  assumptions  about  the  sound  velocity  profile.  Values 
from  20  to  60  seconds  seem  most  reasonable  for  a  range  of  about 
6600  km.  The  duration  scales  directly  with  range. 
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